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TECHNICAL NOTE NO. 581 



A STUDY OF AUTOGIRO ROTOR-BLADE OSCILLA.TIONS 
IN THE PLANE OF THE ROTOR DISK 
By John B. Wheat ley 



SUMMARY 



An analysis of the factors governing the oscillation 
of an autogiro rotor blade in the plane of the rotor disk 
showed that the contribution of the air forces to the re- 
sultant motion was small and that the oscillation is es- 
sentially a direct effect of the rotor-blade flapping mo- 
tion. A. comparison of calculated oscillations with those 
measured in flight on three different rotors disclosed 
that the calculations gave satisfactory agreement with ex- 
periment . The calculated air forces on the rotor blade 
appear to be larger than the experimental' ones , but this 
discrepancy can be attributed to the deficiencies in the 
strip analysis. 



INTRODUCTION 



An autogiro rotor blade is connected to the rotor hub 
in such a way as to permit two articulations, one about 
the horizontal hinge and one about the vertical hinee. 
The horizontal hinge axis is perpendicular to the blade 
span axis and to the rotor axis and permits the blade to 
oscillate freely in a plane containing the blade span axis 
and the rotor axis. The vertical hinge axis is parallel 
to and offset from the rotor axis and permits the blade to 
oscillate in the plane of rotation. This articulation is 
required because the forces acting on the blade in the 
rotor disk are of an unsteady nature and experience has 
shown that heavy stresses and uncomfortable vibrations 
arise if this articulation is not used. In addition, mo- 
tion about this hinge is damped by friction so that tran- 
sient vibrations will quickly subside. 

A knowledge of the motion about the vertical hinge is 
of value for several reasons. The nature of the oscilla- 
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tion must be known in" order to determine intelligently the 
amount of damping required by the. blade. In addition, a 
verified analysis of this oscillation suggests the possi- 
bility of subjecting all rotor characteristics to some 
form of mathematical _a_naly_Bi s and les sens the dependence 
of the des i gner - upon empirical rules-. This paper presents 
an ana-.lysis of this phase of the blade motion and shows a 
comparison between calculated and measured results which 
supports the validity of the analysis. 



ANALYSIS 



The autogiro rotor-blade oscillations in the plane of 
the rotor disk will be divided i'nt-o two parts: Tirst , the 
motion about the vertical pin. that arises from the consid- 
eration of the changing moment of inertia of the bladeB 
about the axis of ro.tat ion caused ..by the flapping motion; 
and second, the forced vibration of the blades about the 
vertical hinge caused by the variation of the air torque 
from a mean value of zero. This method of considering the 
blade motion requires .some, justif ication , because it is 
not obvious that the two oscillations mentioned are inde- 
pendent of each other. 

The experimental " study of vertical-hinge oscillations 
made with a motion-picture camera on the rotor hub estab- 
lished that the amplitude of the oscillation was less than 
0.015 radian and that the osciallation consisted mainly of 
a fundamental of the same frequency as the rotor revolu- 
tions. The angular velocity of thi s . osci 1 lati on is conse- 
quently "less than 0.015 times the mean angular velocity of 
the rotor; therefore the peripheral velocity added by the 
oscillation can safely be neglected, This condition justi 
fies the assumption that the air forces tending to oscil- 
late the blade may be considered independent of this oscil' 
lation. In addition, the amplitude of tho oscillation is 
small enough to justify the assumption that the moment of 
inertia of the blade is not appreciably changed by the os- 
cillation; the effect of the changing moment of inertia 
can, thSr'ef ore, also be considered independent of the ver- 
tical-hinge oscillation. 

Before the analysis Is carried further, it is perti- 
nent to consider briefly the problem of hub vibrations. 
If theblade oscillation about the vertical pin Included a 
harmonic of the same order as the number of blades, the 
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forces arising from the oscillation that acts on the huh 
would add and the motion of the hub would "be affected. 
It can be shown that the sum of the individual effects is 
zero for all harmonics except the one mentioned. It has 
"been experimentally verified that the third and higher or- 
der harmonics are negligible. In addition, experimental 
observations indicate that the hub motion is steady, in- 
asmuch as no steady vibrations in the readings of a sensi- 
tive tachometer were noted. It consequently can be assumed 
that the hub does not oscillate but rotates at a constant 
velocity. . 

Let I be the instantaneous moment of inertia of a 
rotor bla.de and let. I Q be the moment of inertia when the 
angle of flapping p (about the horizontal pin) is zero. 
Then . 

I = I 0 cos 2 p (1) 

The assumption that sin p = j3 is justified by the small 
range covered by p (less than 15°) . Then 

I = I 0 (1 - sin 2 p) = I 0 (1 - p 3 ) (2) 

Let fij be the instantaneous blade angular velocity. The 
law of constant angular momentum states that 

i n ± = Z (3) 



or 

K K 3 

n 4 = - = (1 + p ) (4) 
1 i 0 



where j3 and higher powers are neglected because P is 
small. From reference 1, 



P '= a Q - &i cos — bi sin \j/ - a 2 cos 2 \J/ - b a sin 2 ^ (5) 
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Finally 

- (2a 0 a 1 - a lSs - b^) cos \|/ - (2a Q b x - a x b 2 + b^) sin\)/ 

- (2a Q a s - | ai 8 + | b x 2 ) cos 3* - (2a 0 b a - a l b a ) sin 2t 

- (~a. x a a + b x b a ) cos 3 >J/ - (-a^bg - b iaa ) sin 3 \j/ ] (6) 

Integrating to obtain a displacement equation and letting 
the mean Q i = Q 

0 4 dt = = ~[(1+ a o a + | + +|a s a + iba a )t 

+ n (,3a o b i - a^g + b^^ cos (2a 0 a 1 -a 1 a 2 - b^b a ) sin \J( 

+ i^ (3a o^- a i b ^ °os 8* - (2a Q a 2 - jL-a^ + | b^) sin 3 >|r 

The part of equation (7), varying with ^ which ex- 
presses! the component £ a of the angular displacement of 
the blade oscillation, can be more briefly written as 

t 1 = Uj cos \|f + Vj, sin + | u 3 cos 3 ^ + | 7 8 sin 2 v|/ 

+ 3 u 3 cos 3 ^ + 3 v s sin 3 ^ (8) 

where 

u x = 2 aobi - ai b 2 + b x a a , v x = -2 a^ +■ a 1 a a + b^, 
u 3 = S a Q b a - a 1 b 1 , etc. 
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since 

E 



I 0 Q 1 



The instantaneous air torque acting on a "blade can be 
derived from equation (lO-l) of reference 1 



Q = ^ | P c U 2 cp C L r dr - i P c U 2 8 r dr 



(9) 



The deficiencies in this expression arise from the 
assumptions upon which it is "based. The more important 
ones are: First, that the lift coefficient is proportional 
to the angle of attack and has no maximum value; second, 
that the angle cp, the angle of the local velocity to the 
plane of the rotor disk, is everywhere small enough to 
equate the sine to the angle; third, an average value 8 of 
the "blade profile-drag coefficient is employed which renders 
somewhat uncertain the calculation of instantaneous values 
of the profile-drag torque; fourth., tho lift torque is ex- 
pressed incorrectly where th.8 tangential" velocity Qr is 
less than the forward velocity V; and fifth, the induced 
downward velocity is assumed constant over the disk for 
simplicity, which again renders the derived instantaneous 
values uncertain. The equation is, however, the only avail- 
able means of attacking the problem arid will consequently 
be used, though with caution. 

The solution Q,^ is obtained by substituting for U, 
cp and as in reference 1, and an expression of the 

following type is derived: 

5| = ^ P o a n 2 R 4 (A x cos if + Bi sin * + As cos 2\J/ 

+ B 2 sin 2 + A 3 cos 3 + Bssin 3^ + . . ) (10) 

The coefficients A n and B n are functions of \x, 6, X, 
a n , and b n and are derived in the appendix. The nota- 
tion Q,^ indicates that only the torque varying with \|/ 
is retained, since by hypothesis 
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The component of the angular displacement of the blad 
arising from (fy is designated £ a . If E v is the radi- 
us of the vertical pin, r the radius of the "blade ele- 
ment dr, m x the line density of the blade, R the tip 

radius, and Q E the restoring couple ar is in g fr om £ a , 

then "- - ~ — - 

Q R = / • m x Q 2 R v r £ dr (11) 
s v 

'Che angle i is diagrammed in figure 1. By refer- 
ence to this figure and the law of _sines, assuming small 
angle a , 

• VrX'— 3 - (12) 

■ - ' .... - ' 

Q H = / m x if (r * H T ) R v U dr (is) 

fiy _ , .... , - _ 

R 

It is evident that / m x ( r - R v ) dr integrates 

into the mass moment M m of the "blade about the vertical 
pin - then 

. . Q R = M m n 8 R v ia (14) 

The equation of motion of the blade can now be 
wr itt en 



d t 

I v + M m n a R v U = <ty _ (15) 

where I T is the moment of inertia of the blade about 
the vertical pin. .. . ... ^ — 

No damping term is included because damping of any 



N.A.C.A. Technical Note No. 581 



7 



reasonable magnitude will have a negligible effect on the 
frequency and amplitude of the steady forced -vibration. 

Equation (15) is a well-known form of differential 
equation, the solution of which can be obtained directly 
by assuming 

l s = C 1 cos \|f + D x sin + C 2 cos 2^ + D 2 sin 3* 

+ C 3 cos 3\|/ + D 3 sin 3^ (16) 

The 0 and D coefficients are found by substituting 
for Iq from equation (16) in equation (15) and equating 
the coefficients of trigonometric functions. It is then 
found that 

^ \ p c a R 4 A x ^ 

D = 1 P ° & R4 Bl (is) 

£ p c a R 4 A a 

C a = (19) 

"m ~~ ^o 

i p c a R 4 B a 

p c a R 4 A 3 . 
°3 = M m R v ~- 9I 0 Ul) 

p c a R 4 B, 

^ = 1_ 3 (22) 

Mm R y — 9 I 0 

The resultant blade displacement £ can now be ob- 
tained by adding equation (16) to equation (8), neglecting 
all constant terms, then 

t, = Ej_ cos \j/ + I" x sin + E a cos 2\!/ + F a sin 2\(/ 

+ E 3 cos 3^ +' F 3 sin 3^ (23) 

where 
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Pi = D x + Ul 
TS a = C 3 + | v a 
*a = D a + | u 2 



(24) 



= C, 



*3 = 3 3 + 



3 
1 



u. 



It is of interest to consider a: simplified solution 
in which the C and D eg ef f^ici^ent s a.r.e.Jn_eglected and a 
flapping curve is assumed that has no second harmonics. 
Then 



p = a 0 - a x cos \|/ - \> 1 sin i(r 
a 0 +. -v/a^ T \ 3 



max 



(25) 
(26) 



£ *3 cos \(/ - Sa^ sin ij. a x bj cos 3 \J/ 



Tho maximum value of £ is, noglecting the second 
harmonics, which are normally quite small: 

^max = 2 a o ^ + V = 2 * 0 (iW - * 0 > 



(38) 



TESTS AND APPAEA TUS 



The rotor-blade motions about the horizontal and ver- 
tical pins of several autogiros were measured by mounting 
a motion-picture camera on the rotor hub and photographing 
the rotor blade in flight (reference 3). Data simultaneous 
ly obtained determined the tip-speed ratios. One rotor, 
the PCA-3, was also tested alone in tho full-scale wind tun 
nel (reference 3), the tests affording somo necessary data 
shown in_ subsequent tables. The physical characteristics 
of the three rotors tested, the PCA-3, the PAA-1 , and the 
KD-1, are given" in tables I, II, and III, respectively. 



RESULTS 



The flapping coefficients of the PCA-3 rotor are pre- 
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sented as functions of the tip-speed ratio in figure 2. 
Figures 3, 4, and 5 present comparisons of the experimental 
vertical-pin oscillation in coefficient form with values 
calculated from the air forces alone, the angular momentum 
alone, and the combined air forces and angular momentum, 
respectively. No third harmonics are givon since both cal- 
culated and experimental values were found to he negligible. 
A comparison of the vertical-pin motion during a rotor rev- 
olution at a tip-speed ratio of 0 ..5 , based on calculated 
and experimental results, is shown in figure 6. 

The experimental flapping coefficients of the PAA-1 
rotor are presented in figure 7 as functions of the tip- 
speed ratio. The coefficients of the vertical-pin oscilla- 
tion of this rotor based on experiment and upon calculation 
from the angular-momentum law are given in figure 8. 

Figure 9 shows the experimental flapping coefficients 
of the KD-1 rotor as a function of tip-speed ratio. The ex- 
perimental coefficients of the vertical-pin oscillation of 
this rotor are compared in. figure lo with the coefficients 
calculated from the angular momentum alone. 

All calculations are based upon the data in tables I, 
II, and III and in figures 2, 7, and 9, which show the ex- 
perimental blade flapping-motion data for each rotor. 

DISCUSSION 

In order that the desiener be able to calculate the 
rotor-blade motion about the vertical pin, he must first 
calculate the flapping motion about the horizontal pin. 
This problem has been analyzed and discussed in reference 
1 and will not be touched upon here. .. 

The most important deduction in this study is arrived 
at from an examination of figures 3, 4, and 5, which demon- 
strate that the blade motion about the vertical pin arising 
from the air forces is negligible compared with the motion 
caused by the angular momentum. In addition, figures 4, 8, 
and 10 indicate that a reasonably accurate calculation of 
the blade motion can be made without considering the air 
forces . 

The agreement between the calculated and experimental 
blade-motion coefficients is individually excellent as shown 
in figures 4, 8, and 10, except for the component u x , the 
sin \j/ term. The poor agreement between this calculated 
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factor and the measured one can "be explained by figure 3_ t 
The factor £>i represents the calculated sin component 
arising from the air forces. . It appears that the influence 
of the air forces on the motion is reflected, almost entire- 
ly in the coefficient of sin \J/ and that the D x term rep- 
resenting this effect is of the proper sign to make m 1 
agree more closely with experiment. It is found, however, 
in figure 5. that , which is the sum of u x and D x , dis- 
agrees with, the experiment as radically as does v. x . It 
must therefore he concluded that the calculation has over- 
estimated the magnitude of D x . This deduction is reason- 
able be" caus e the actual accelerating torque on the retreat- 
ing side of the disk must be less than the calculated torque 
bocause the strip analysis assumes a constant lift-curve 
slope and does not take into account the existence of a max- 
imum lift; both of these factors wsuld cause an ov.er-estima- 
tion of the accelerating torque where the local angles of 
attack are large. 

The plotted blade motion in figure 6 illustrates graph- 
ically the error in the motion arising from the u x term; if 
this calculated term agreed with the experiment the two 
curves of figure 6 would, for all practical purposes, coin- 
cide. ' ' ' 

An examination of figures 4, 8, and 10 suggests that a 
very cLo-Sjs approximation to the. actual motion about the ver- 
tical pin can be obtained by calculating the motion from the 
angular momentum alone and arbitrarily addine to the Uj, 
term a correction factor of the magnitude indicated in these 
f i eur.es . 

- ' CONCLUSIONS 

1. The autogiro rotor-blade oscillation about the ver- 
tical pin is essentially a direct effect of the rotor-blade 
flapping motion. 

2. The effect of the air forces on the autogiro rotor- 
blade oscillation about the vertical pin was shown by both 
calculation and experiment to be of second order in compari- 
son to the effect of the blade flapping. 

3. The deficiencies of the strip analysis result in the 
overestimation of the air forces acting on the rotor blades. 



Langley Memorial Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va . , June 23, 1936. 
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APPENDIX 



from reference 1, equation (lO-l). 

BR R 
Q = y ^pcU 2 qpO L rdr- / ^ p c U a 6 r dr 

o " o 
Substitute for U, cp and C L as in reference 1; 

integrate, then, dropping all constant terms that must sum 
to zero since the mean rotor angular velocity is constant, 

^ = |pcaQ 3 R 4 {(^a 0 [ |e o B 3 + i 6l B 4 ] + b x [ |e o B 4 + l 01 B 5 _ 

-^a 0 B a +|Xh x B 3 +1^X32 B 2 - |ia Q ba B 3 
-|j2a 0 a 1 B 2 +|iJLa 1 b 1 B 3 +|a 1 a s B 4 +|b 1 b s B 4 ) cos ^ 

, + ( l xxr|e 0 B 2 + |e 1 B 3 ]-a 1 TleoB-H-le.B 5 ] + ^ 3 ai [^b^b 3 ] 

- ^ba [|9 o B 3 +| 6x B 4 ] - |^a x B 3 +|m> b a B a - Jp 8 a 0 b x B 8 
+ l^aoaa B 3 - ip*^ B 3 +|^ 1 ba B^b! 8 B 3 - |b, a s B 4 - sin * 

• 4(p.a x [|e o B 3 + | 01 B 4 ] +2b a [leoB 4 -.^, B 5 ] a, B 8 

+|xba B 3 +^ a a 0 S B 8 - ^a^ B 3 -§y, 3 a^ B 3 -^ 3 B 4 
+ |y 3 a 1 3 ^+|ia 1 b 3 B 3 +|-b 1 3 B^lnbiasB 3 ^-^ 3 !) cos 2* 
+(-JpRa 0 [leoB 2 ^^ 3 ]^ [ie 0 B 3 + |6 1 B 4 J-2a a \^ Q B* + U^ 
+ |y.X bi B 3 - 1?\ a 2 B 3 +|y,a 0 ax B 3 - gp, 8 a e b a B 3 
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- laibiB 4 +1^^x01 B 2 -|ua ia2 B 3 bibs B 3 ^ Bin 2* 
-^(^a^^^b.rieoB^le.B 3 ] + |aa 2 [|e 0 B 3 + | 9l B*j 

- |u 3 a 0&1 B a -|txa 0 b3 B 3 +| i aa 1 b a B 3 a 2 B 4 4™bi ha B 4 ) cos 3\{/ 

+ (i^aB^^a^^^ 

+ B 3 -|b 1 a s B 4 ^ sin 
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TABLE I. PGA-2 Rotor Characteristics 



Radius E, ft 22.5 

Moment of inertia about vertical pin I x , slug-ft 316 

Radius of vertical pin R v , ft .646 

Blade mass m, slugs 2.588 

Blade chord (outer part of blade) c, ft 1.833 

Lift-curve Blope a 5.85 

Blade span efficiency factor B .959 

Blade pitch angle 6, radians -0.0157 + 0.0872 ~ 



(The values of ^ and 6/a are calculated from reference 3) 
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TABLE II. PAA-1 Rotor Characteristics 

Radius R, ft 19 

Estimated moment of inertia about 

vertical pin I x . slug-ft? 150 

Estimated radius of vertical pin R v , ft. 0.75 

Blade mass m, slugs 1.65 

TABLE III. ZD-1 Rotor Characteristics 

Radius A, ft 20 

Estimated moment of inertia ahout 

vertical pin I x , slug-ft. s 155 

Estimated radius of vertical pin R v , ft 0.792 

Blade mass m, slugs ' . 1.61 




Figure 1.- Geometry of autogiro rotor "blade and vertical-pin 
articulation. 
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Figure 3.- Coefficients of PCA-§ rotor-Made motion about the vertical 

pin as measured in flight and as calculated from the air 
forces alone. 
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Figure 4.- Coefficients of PCA.-2 rotor-"blade motion about the vertical 

pin as- measured in flight and as calculated from- the angular 
momentum alone . 
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Figure 5.- Coefficients of PCA-2 rotor-blade motion about the vertical 

pin as measured in flight and as calculated from the air 
forces and the angular momentum. 
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Figure 3.- Coefficients of the PAA-1 rotor-olade motion about the 

vertical pin as measured in flight and as- calculated 
from the angular momentum alone. 




Figure 9.- 



Flapping coefficients of the KD-1 autogiro rotor blade 
as measured in flight. ~ 
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Figure 10.- Coefficients of the KD-1 rotor-hlade motion ahout the 

vertical pin as measured in flight and as calculated 
from the angular momentum alone. 



